A Search for Periodicities in the Solar Neutrino Flux 
Measured by the Sudbury Neutrino Observatory 
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A search has been made for sinusoidal periodic variations in the *B solar neutrino flux using data 
collected by the Sudbury Neutrino Observatory over a 4-year time interval. The variation at a period 
of one year is consistent with modulation of the *B neutrino flux by the Earth's orbital eccentricity. 
No significant sinusoidal periodicities are found with periods between 1 day and 10 years with either 
an unbinned maximum likelihood analysis or a Lomb-Scargle periodogram analysis. The data are 
inconsistent with the hypothesis that the results of the recent analysis by Sturrock et al, based on 
elastic scattering events in Super-Kamiokande, can be attributed to a 7% sinusoidal modulation of 
the total *B neutrino flux. 

PACS numbers: 26.65+t, 95.75.Wx, 14.60.St, 96.60.Vg 



I. INTRODUCTION 



There have been recent reports of periodic variations 
in the measured solar neutrino fluxes [1-8]. Other anal- 
yses of these same data, including analyses by the ex- 



perimental collaborations themselves, have failed to find 
such evidence [9, 10]. The reported periods have been 
claimed to be related to the solar rotational period. Par- 
ticularly relevant for this paper is a claimed 7% ampli- 
tude modulation in Supcr-Kamiokandc's ^B neutrino flux 
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at a frequency of 9.43 y""'^ [2, 3]. Because solar rotation 
should not produce variations in the solar nuclear fusion 
rate, non-standard neutrino properties have been pro- 
posed as an explanation. For example, the coupling of 
a neutrino magnetic moment to rotating magnetic fields 
inside the Sun might cause solar neutrinos to transform 
into other fiavors through a resonant spin fiavor preces- 
sion mechanism [11-13]. Periodicities in the solar neu- 
trino flux, if confirmed, could provide evidence for new 
neutrino physics beyond the commonly accepted picture 
of matter-enhanced oscillation of massive neutrinos. 

This paper presents a search for periodicities in the 
data from the Sudbury Neutrino Observatory (SNO). 
SNO is a real-time, water Cherenkov detector located in 
the Inco, Ltd. Creighton nickel mine near Sudbury, On- 
tario, Canada [14]. SNO observes charged-currcnt (CC) 
and neutral-current (NC) interactions of neutrinos on 
deuterons in 1 ktonne of D2O, as well as neutrino-electron 
elastic scattering (ES) interactions. By comparing the 
observed rates of CC, NC, and ES interactions, SNO has 
demonstrated that a substantial fraction of *B electron 
neutrinos produced inside the Sun transform into other 
active neutrino flavors [15-19]. 

SNO's combination of real-time detection, low back- 
grounds, and sensitivity to different neutrino flavors give 
it unique capabilities in a search for neutrino flux period- 
icities. Chief among these is the ability to do an unbinned 
analysis, in which the event times of individual neutrino 
events are used as inputs to a maximum likelihood fit. 

This paper presents results from an unbinned maxi- 
mum likelihood analysis and a more traditional Lomb- 
Scargle periodogram analysis for SNO's pure D2O and 
salt phase data sets. Previous analyses of data from 
other experiments have used the Lomb-Scargle peri- 
odogram [20] and binned maximum likelihood techniques 
to search for periodicities in the solar neutrino data. 
These data generally consist of flux values measured in 
a number of time bins of unequal size. Because analyses 
of binned data can be sensitive to the choice of binning, 
which can also produce aliasing effects, it is desirable to 
avoid binning the data if possible. Section II describes 
the data sets. Section III contains the results of a gen- 
eral search for any periodicities with periods between 1 
day and 10 years. Section IV presents limits on the am- 
plitudes at two specific frequencies: the 9.43 y~^ mod- 
ulation of the *B neutrino flux claimed by Sturrock et 
al. [2, 3], and a yearly modulation due to the Earth's 
orbital eccentricity. 



II. DESCRIPTION OF THE SNO DATA SETS 

The data included in these analyses consist of the se- 
lected neutrino events for the initial phase of SNO, in 
which the detector contained pure D2O [16], and for 
SNO's salt phase, in which 2 tonnes of NaCl were added 
to the D2O to increase the neutron detection efficiency for 
the NC reaction [19]. Each data set is divided into runs 



of varying length during which the detector was live for 
solar neutrino events. The D2O data set consists of 559 
runs starting on November 2, 1999, and spans a calen- 
dar period of 572.2 days during which the total neutrino 
livetime was 312.9 days. The salt phase of SNO started 
on July 26, 2001, 59.7 calendar days after the end of the 
pure D2O phase of the experiment. The salt data set con- 
tains 1212 runs and spans a calendar period of 762.7 days 
during which the total neutrino livetime was 398.6 days. 
The intervals between runs during which SNO was not 
recording solar neutrino events correspond to run tran- 
sitions, detector maintenance, calibration activities, pe- 
riods when the detector was off, etc. Deadtime incurred 
within a run, mostly due to spallation cuts that remove 
events occurring within 20 seconds after a muon, can be 
neglected, since such deadtime is incurred randomly at 
average intervals much shorter than the periods of inter- 
est for this analysis. This deadtime is 2.1% for the D2O 
data set and 1.8% for the salt data set. 

The event selection for the data sets is similar to that in 
[16] and [19]. Events were selected inside a reconstructed 
fiducial volume of i? < 550 cm and above an effective 
kinetic energy of T^s > 5 MeV (D2O) or Toff > 5.5 MeV 
(salt). The salt data set contains 4722 events, as in [19]. 
During the salt analysis described in [19] a background 
of "event bursts" , consisting of two or three neutron-like 
events occurring in a short time interval, was identified 
and removed with a cut that eliminated any event occur- 
ring within 50 ms of an otherwise acceptable candidate 
neutrino event. The source of these 11 bursts is not cer- 
tain, but they may have been produced by atmospheric 
neutrino interactions. For this analysis a similar cut re- 
moving any event occurring within 150 ms of another 
event was applied to the D2O data, reducing the number 
of selected events from 2928, as in [16], to 2924. The tim- 
ing window for the cut in the D2O data is longer than for 
the salt data to account for the longer neutron capture 
time in pure D2O. 

An important element of a periodicity analysis is ex- 
act knowledge of when each data-taking run began and 
ended. These run boundaries define the time exposure 
of the data set, which itself may induce frequency com- 
ponents that could impact a periodicity analysis. The 
unbinned maximum likelihood analysis described below 
makes explicit use of these run boundary times, and all 
Monte Carlo simulations arc generated using the exact 
run boundaries, even if the simulated data are binned 
in a following analysis. These precautions avoid ad hoc 
assumptions about the distribution of the time exposure 
within any time bin. The measured time for each event 
was measured with a global positioning system (GPS) 
clock to a precision of ~ 100 ns, but rounded to 10 ms 
accuracy for the analysis. The run boundary times were 
determined from the times of the first and last events in 
each run with a precision of ~ 50 ms. 

Figure 1 displays the solar neutrino event rate in live- 
time corrected 1-day bins over the total exposure time 
of both phases of SNO [21]. The D2O and salt data sets 
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FIG. 1: The livetime-corrected 1-day total rate of events as a 
function of time for the D2O and salt phases. The weighted 
mean rates are 9.35 ± 0.17 and 11.85 ± 0.17 events/day, re- 
spectively. 



may be individually examined for periodicities, or the 
combined data from both phases can be jointly searched. 
It should be noted that the relative amounts of CC, NC, 
and ES events are different for the D2O and salt data, 
with the salt data set containing a much higher fraction 
of NC events. 

Although SNO's data sets are dominated by solar 
neutrino events, they also contain a small number of 
non-neutrino backgrounds, primarily neutrons produced 
through photodisintegration of deuterons by internal 
or external radioactivity. The total estimated number 
of background events is 1231^7 ^^'^ D2O data set 
{4.2to,l% of the total rate), and 260±59 events for the 
salt data set (5.5±1.2% of the total rate). Although the 
background rate is not entirely constant, the backgrounds 
are small and stable enough that they can be neglected 
in this analysis. 



III. GENERAL PERIODICITY SEARCH 

Both an unbinned maximum likelihood analysis and a 
Lomb-Scarglc pcriodogram with 1-day binning were used 
to search SNO's data for periodicities. Results arc pre- 
sented below for the D2O, salt, and combined data from 
each method, along with evaluations of the sensitivity of 
each method to sinusoidal variations of various periods 
and amplitudes. The periodicity searches were carried 
out over the sum of CC, NC, and ES events. 

Extensive use was made of Monte Carlo data sets to 
evaluate the statistical significance of the results and the 
sensitivity of each method. To determine the statisti- 
cal significance of any peak in the frequency spectrum, 
10,000 Monte Carlo data sets with events generated ran- 
domly within the run boundaries for each phase were 
used, with mean event rates in each phase matching those 



observed in SNO's data sets. The number of events in 
each Monte Carlo data set was drawn from a Poisson dis- 
tribution with the same average rate as the data, and the 
events were distributed uniformly within the run bound- 
aries [22]. These "null- hypothesis" Monte Carlo data 
sets were used to determine the probability that a data 
set drawn from a constant rate distribution would pro- 
duce a false positive detection of a periodicity. To deter- 
mine the sensitivity of an analysis to a real periodicity, 
1,000 Monte Carlo data sets were generated for each of 
several combinations of frequencies and amplitudes, with 
the events drawn from a time distribution of the form 
1 -I- Acos(2Tr ft). The sensitivity for any frequency / and 
amplitude A is then defined as the probability that the 
analysis will reject the null hypothesis of a constant rate 
at the 99% confidence level. 



A. Unbinned Maximum Likelihood Method 

The unbinned maximum likelihood method tests the 
hypothesis that the observed events are drawn from a 
rate distribution given by 

(t){t)=N{l + Acos{2nft + S)) (1) 

relative to the hypothesis that they are drawn from a 
constant rate distribution (A = 0). A is the fractional 
amplitude of the periodic variation about the mean, S 
is a phase offset, and is a normalization constant for 
the rate. Equation 1 serves as the probability density 
function (PDF) for the observed event times, which are 
additionally constrained to occur only within run bound- 
aries (i.e., (j){t) = if t is not between the start and end 
times of any run). 

With / fixed, the maximum of the extended likelihood 
L{tk\N, A, S, f) as a function of the individual event times 
tk is calculated for a data set as 

runs events 

\nL{tk\N,A,SJ) = -J2 'mdt+ ^ ln(0(tfc)) 
3=1 k=i 

(2) 

where the first term is a sum over all runs of an inte- 
gral evaluated between each run's start and stop times 
tl and and accounts for Poisson fluctuations in the 
signal amplitude. The second term is a sum over the 
events in the data set, and tk is the time of the fcth 
event. The log likelihood is maximized as a function of 
A, S, and N to yield Inimax, while / is kept fixed. Then 
the constraint A = is imposed, removing the depen- 
dence of (f){t) on both A and 6, and the log likelihood is 
maximized over the remaining free parameter N to yield 
lnLniax(^ = 0). By the likelihood ratio theorem [23] the 
difference 25 = 2(lniniax^liiLmax(^ = 0)) will approxi- 
mately have a distribution with two degrees of freedom 
(since the choice A ~ Q also removes the dependence on 
the phase 6). Thus 5* will follow a simple exponential 
if the true value of A is zero. Therefore, at any single 
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frequency /, the probability of observing S > Z under 
the null hypothesis that 4){t) ~ N \s approximately . 
This null hypothesis test is carried out for a large set of 
frequencies scanning the region of interest. 

Equation 2 includes both a floating offset N and an 
amplitude A as free parameters. Allowing both of these 
parameters to vary is necessary to deal with very low 
frequencies, for which N and A become degenerate pa- 
rameters. Simply fixing TV to the mean rate, as was done 
in [3], will be prone to bias at the very lowest frequencies, 
but gives virtually identical results to the floating offset 
procedure when the length of the data set is longer than 
the period T = 1//, since in this case enough cycles are 
sampled to break the degeneracy between A'' and A. 

Equations 1 and 2 are adequate to test for periodicity 
in a single data set, but for a combined analysis of SNO's 
D2O and salt data sets, account must be taken of the dif- 
fering mean rates owing to different detection cfRciencies 
and energy thresholds in the two phases. This can be 
done by generalizing (f>{t) to: 

(i>{t) = iVd2o(l + ^cos(27r/i + (5)), lite D2O run 
(t>{t) = iVsait(l + Acos(27r/t + 5)), if t e salt run 

This PDF allows different normalization constants for the 
two data sets, while retaining the assumption that the 
flux variation has the same fractional amplitude in both 
the D2O and salt data. 
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FIG. 2: Log likelihood difference [S) as a function of fre- 
quency for the unbinned maximum likelihood method for 
SNO's D2O, salt, and combined data sets. 



1. Results for the SNO data sets 

Figure 2 shows 5 as a function of frequency for the 
D2O, salt, and combined data sets at 3650 frequencies 
with periods ranging from 10 years down to 1 day, with 
a sampling interval of A/ = 1/(3650 days). This corre- 
sponds to an oversampling of the number of independent 
Fourier frequencies for continuous data by a factor of ap- 
proximately 5-6 for the separate D2O and salt data sets, 
and a factor of 2.6 for the combined data set. The max- 
imum value of S for the D2O data set is 5 = 7.1 at 
a period of 3.50 days (/ = 0.296 days^^). The largest 
peak found in the salt data has a height of S* = 9.4 at a 
period of 1.03 days (/ = 0.971 days~^), while the com- 
bined data set has its largest peak of S* = 8.8 at a period 
of 2.40 days (/ 0.417 days'i). 

Figure 3 shows the distribution of maximum peak 
heights for 10,000 Monte Carlo data sets generated with 
no periodicity and analyzed identically to SNO's com- 
bined data set. Distributions for the D2O and salt Monte 
Carlo data sets analyzed individually look similar. Of the 
10,000 simulated data sets, 35% yielded at least one peak 
with S > 8.8, exceeding the largest peak seen in SNO's 
combined data set. For the D2O Monte Carlo data sets, 
72% had a peak larger than the observed largest peak 
of S' = 7.1, while 14% of the salt Monte Carlo data sets 
yielded a peak larger than the S = 9.4 peak seen in the 
data. Therefore, none of the observed peaks are statisti- 
cally signiflcant. 
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FIG. 3: Distribution of the maximum value of S for 10,000 
Monte Carlo data sets produced with A = for a combined 
D20H-salt unbinned maximum likelihood analysis. A fit of 
the distribution to Equation 3 is shown. 
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Under the null hypothesis of no time variability, the 
probability of any individual frequency having S smaller 
than some threshold Z is approximately 1 — e^^. If all 
3650 scanned frequencies were statistically independent, 
the probability that all peaks would be smaller than Z 
would be (1 — e^^)"^^^". However, the 3650 scanned fre- 
quencies arc not strictly independent, since a finite data 
set has limited frequency resolution, and neighboring fre- 
quencies are correlated. If F is the effective number of 
independent frequencies, then the probability distribu- 
tion for the height Z of the largest peak approximately 
follows 



P{Z) dZ oc e-^il 



-Z\F-1 



dZ 



(3) 



The effective number of independent frequencies in- 
creases with the length of the data set and number of 
detected events. Fitting the Monte Carlo distributions 
for Z to this equation yields F = 1422 ± 17 for the 
D2O, F = 1696 ± 26 for the salt, and F = 2739 ± 30 
for the combined data set. Figure 3 shows this fit for the 
combined analysis. These values are consistent with ex- 
pectations based on the oversampling factors described 
in Section III A 1 [24] . Although Equation 3 appears to 
model P{Z) well, quoted significance levels are always 
determined directly from the Monte Carlo distributions 
and not from the analytic formula. To ensure that no 
significant peaks were missed, the combined analysis of 
the actual data (but not the Monte Carlo data sets) was 
repeated with the sampling increased by a factor of five. 
No new peaks were found. 



2. Sensitivity to sinusoidal periodicities 

Distributions of the maximum peak height for Monte 
Carlo data sets, such as in Figure 3, readily yield the 
threshold C for which 99% of Monte Carlo data sets 
generated without periodicity would yield a maximum 
peak height of C, or less. This threshold defines the peak 
height at which the null hypothesis of no time variation 
is rejected at the 99% confidence level, and equals 12.10, 
12.20, and 12.65 for the D2O, salt, and combined data 
sets respectively. 

Monte Carlo data sets drawn from rate distributions 
with sinusoidal periodicities of various periods and am- 
plitudes were analyzed to determine the probability of 
rejecting the null hypothesis at the 99% C.L. 

Figure 4 shows the amplitudes as a function of period 



at which the method has a 50% (90%) probability of re- 
jecting the null hypothesis, for simulations of SNO's com- 
bined data set. While the sensitivity varies as a function 
of period, a signal must have an amplitude of approxi- 
mately 8% to be discovered 50% of the time. 
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FIG. 4: Sensitivity contours for the unbinned maximum like- 
lihood and Lomb-Scargle methods. The contours indicate 
the amplitude as a function of period at which the analysis 
method has a 50% (90%) chance of discovering a sinusoidal 
variation of unknown period at the 99% confidence level. For 
longer periods both methods are equally sensitive. Due to its 
lack of binning, the maximum likelihood method has sensitiv- 
ity to much smaller periods than the Lomb-Scargle method. 



B. The Lomb-Scargle periodogram 

The Lomb-Scargle periodogram is a method for search- 
ing unevenly sampled data for sinusoidal periodicities [20] 
and provides an alternative to the unbinned maximum 
likelihood technique described above. 

The Lomb-Scargle power P{f ) at frequency / is calcu- 
lated from the measured flux values y{ti) in N indepen- 
dent time bins as: 



J 



P{I) 




E^i WtiviU) - y) cos(27r/(t, - t)) 



y) sin(27r/(t, - r)) 



Eili cos2(27r/(t, - r)) 



(4) 
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where the phase factor t satisfies: 



tan(47r/7 



Eili ■'^i cos(47r/tj) 



Each bin is weighted in proportion to the inverse of 
its squared uncertainty divided by the average value of 
the inverse of the squared uncertainty (l/crf) (so Wi = 
1/af / (l/af)), as in [3]. In Equation 4 is the livetime- 
weighted mean time for the ith bin, and y and are the 
weighted mean and weighted variance of the data for all 
the bins calculated with the weighting factors Wi. 

Like the maximum likelihood method, the power P in 
the Lomb-Scargle periodogram at any single frequency / 
is expected to approximately follow an exponential dis- 
tribution if the data set is drawn from a constant 
rate distribution. The same methods of evaluating the 
significance of the largest peak and the sensitivity of the 
method to periodic signals can be employed, making use 
of large numbers of Monte Carlo data sets. 

In [9] the Supcr-Kamiokandc collaboration used an un- 
weighted Lomb-Scargle periodogram {wi = 1) to search 
its data set for periodicities, a choice that was criticized 
in [3]. The analysis presented here used the weighted 
Lomb-Scargle periodogram. 

For the Lomb-Scargle method SNO's recorded events 
were binned in 1-day intervals (see Figure 1), and the 
livetime, the livetime- weighted mean time ti, and the 
event rate y{ti) were calculated for each bin. To pre- 
vent biases stemming from the assumption of Gaussian 
statistics, any bin in which fewer than five events would 
be expected based upon that bin's livetime and the mean 
event rate was combined with the following bin(s) so that 
the expected number of events in all bins was greater than 
five. The uncertainty ai on the rate in each bin was taken 
to be the square root of the expected number of events 
in that bin for a constant rate. This calculation of the 
uncertainty is appropriate if one views the Lomb-Scargle 
method as a null hypothesis test of the no-periodicity hy- 
pothesis; however, using the observed number of events 
instead to calculate ai does not change the conclusions 
of this study. 

When doing a combined D2O -I- salt analysis one must 
account for the different mean event rates in the two 
phases. In the Lomb-Scargle analysis this was accom- 
plished by scaling the rates and uncertainties on the salt 
data bins by the ratio of the weighted mean D2O rate to 
the weighted mean salt rate. 
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FIG. 5; Lomb-Scargle periodograms for SNO's D2O, salt, and 
combined data sets, with 1-day binning. 



cies less than 0.5 days ^ to avoid potential binning ef- 
fects. The maximum peak height for the D2O data set 
is P = 7.1 at a period of 2.45 days (/ = 0.408 days"^). 
The largest salt peak has a height of P = 6.8 at a pe- 
riod of 2.33 days (/ = 0.429 days^^), while the combined 
data set has its largest peak of P = 8.7 at a period of 
2.42 days (/ = 0.413 days'i). 

The probability of observing a larger peak than that 
actually seen in the Lomb-Scargle periodogram, if the 
rate were constant, was estimated using the previously 
described 10,000 Monte Carlo data sets having no peri- 
odicity. Under the null hypothesis of no time variability, 
the probability of getting a peak larger than the biggest 
peak seen in the Lomb-Scargle periodogram is 46% for 
the D2O data set, 65% for the salt data set, and 27% 
for the combined data sets. As with the unbinned maxi- 
mum likelihood method, no evidence for time variability 
is seen. 



2. Sensitivity to sinusoidal periodicities 



1. Results for the SNO data sets 

Figure 5 shows the Lomb-Scargle periodograms for 
the D2O, salt, and combined data sets. A total of 
7300 frequencies were tested with periods ranging from 
10 years to 2 days, with a samphng step of A/ = 
1/(14600 days) [25]. Because the data were binned in 
one-day intervals, the analysis was restricted to frequen- 



Monte Carlo data sets generated with sinusoidal pe- 
riodicities were used to estimate the sensitivity of the 
Lomb-Scargle method to signals of various periods and 
amplitudes. Figure 4 shows the amplitudes as a func- 
tion of frequency at which the analysis would detect the 
signal 50% or 90% of the time. In each case the signal 
is considered to be detected if the Lomb-Scargle method 
rejects the null hypothesis of a constant rate at the 99% 
confidence level. The threshold for rejecting the null hy- 
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FIG. 7: Relative event rate as a function of days since perihe- 
lion, normalized to the mean rate. In this plot the mean D2O 
event rate has been scaled to match the mean salt phase event 
rate, and the two data sets have been combined. The curve 
represents the expected variation due to the eccentricity of 
the Earth's orbit. 



FIG. 6: Distribution of maximum Lomb-Scargle powers for 
1,000 Monte Carlo data sets produced with amplitudes A — 
0.0, 0.10, 0.15, 0.20, 0.25 and a period of 20 days for the com- 
bined D2 0-|-salt analysis. 



pothesis at the 99% C.L. is P = 11.15 for the D2O data, 
P = 11.43 for the salt data, and P = 12.24 for the com- 
bined analysis. Figure 6 shows example maximum power 
distributions for a 20-day period with amplitudes of 0, 
10, 15, 20, and 25% for the combined analysis. 

3. Systematic checks 

Many checks of the Lomb-Scargle periodogram were 
made to verify that the results are robust. In particular, 
all data and Monte Carlo results were recomputed for (a) 
a range of bin sizes, from 1-day to 5-days in fractional 
day steps, (b) a range of starting times of the first bin 
in fractional day steps, and (c) different values of the 
frequency sampling step. There was no evidence for time 
variability under any of these scenarios. 

IV. LIMITS AT SPECIFIC FREQUENCIES OF 
INTEREST 

The sensitivity calculations in Sections IIIA2 and 
III B 2 are appropriate when the frequency of the sig- 
nal is not known a priori^ and could occur anywhere in 
the frequency search band. The threshold for claiming a 
detection at the 99% C.L. must accordingly be set rela- 
tively high to reduce the false alarm probability, which 
was found from Monte Carlo simulations but is approx- 
imately given by integrating Equation 3 above the de- 
tection threshold, to < 1%. However, if the frequency 
of interest is specified a priori^ then a more restrictive 



and sensitive test can be done using the fitted amplitude 
at that frequency. Two particular frequencies of interest 
are the 7% variation in the Super-Kamiokande data at 
a frequency of 9.43 y^^ claimed by Sturrock et al. [3], 
and the annual modulation of the neutrino flux by the 
Earth's orbital eccentricity. 

A. Test at / = 9.43 y"^ 

Sturrock et al. have claimed evidence for a periodic- 
ity in Super-Kamiokande's neutrino data at a frequency 
of 9.43 ± 0.05 y-i (0.0258 days'^) with an amplitude of 
7% [3]. Examination of SNO's unbinned maximum like- 
lihood results in the interval from 9.33-9.53 y""'^ yielded 
no value larger than S = 1.1 in either the D2O, the salt, 
or the combined data sets [26]. The best-fit amplitude 
for the combined data set inside this frequency interval 
is (1.3 ± 1.6)%. This disagrees with a 7% amplitude pe- 
riodicity in the neutrino flux by 3.6 sigma. It must 
be remarked that SNO's limit applies to a modulation 
of the summed rates of CC, ES, and NC events above 
their respective energy thresholds, whereas the reported 
7% periodicity in the Super-Kamiokande data is a mod- 
ulation of the elastic scattering rate from ®B neutrinos 
above a total electron energy threshold of 5 MeV. The 
best- fit amplitudes for the D2O and salt data sets are 
(3.8 ± 2.6)% and (0.3 ± 2.3)% respectively. 

B. Eccentricity Result 

The Earth's orbital eccentricity is expected to produce 
a rate variation proportional, in excellent approximation, 
to (1 ecos(27r/(t - ^o)))^ where e = 0.0167 is the ec- 
centricity of the orbit, / = 1/(365.25 days) is the Earth's 
orbital frequency, and tg is the time of perihelion. Max- 
imum sensitivity to this effect is obtained if to and / are 
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fixed to their known values and the combined data sets 
are fit for e only. This has been implemented using the 
unbinned maximum likelihood technique. The best-fit 
eccentricity is e = 0.0143 ± 0.0086, in good agreement 
with the expected value. The difference in the log like- 
lihoods for the best fit compared to e = is 1.394. The 
probability of obtaining a larger value of the log likeli- 
hood difference if e = is 9.5%. Figure 7 displays the 
relative event rate for the combined data as a function of 
the time since perihelion. 

V. CONCLUSIONS 

Data from SNO's D2O and salt phases have been exam- 
ined for time periodicities using an unbinned maximum 
likelihood method and the Lomb-Scargle periodogram. 
No evidence for any sinusoidal variation is seen in either 
data set or in a combined analysis of the two data sets. 
This general search for sinusoidal variations with periods 
between 1 day and 10 years has significant sensitivity 
to periodicities with amplitudes larger than ^ 8%. The 
best-fit amplitude for a sinusoidal variation in the total 
neutrino flux at a frequency of 9.43 y~^ is (1.3±1.6)%, 
which is inconsistent with the hypothesis that the results 
of the recent analysis by Sturrock et al. [3], based on 



elastic scattering events in Super-Kamiokande, can be at- 
tributed to a 7% modulation of the neutrino flux. A 
fit for the eccentricity of the Earth's orbit from the mod- 
ulation at a period of one year yields e = 0.0143 ±0.0086, 
in good agreement with the known value of 0.0167. 
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